arXiv:astro-ph/0512160v2 29 Aug 2006 


Mon. Not. R. Astron. Soc. 000 , 000-000 (0000) Printed 5 February 2008 (MN style file v2.2) 


Cluster infall in the concordance A-CDM model. 


Maximiliano C. Pivato 1 , Nelson D. Padilla 1,2 and Diego G. Lambas 1,3 

1 Grupo de Investigaciones en Astronomia Teorica y Experimental (IATE), Observatorio Astronomico Cordoba, UNC, Argentina. 

2 Departamento de Astronomia y Astrofisica, Pontificia Universidad Catolica, V. Mackenna 4860, Santiago 22, Chile. 

3 Consejo Nacional de Investigaciones Cientificas y Tecnologicas (CONICET), Argentina. 


Accepted ???? Month ??. Received ???? Month ?? 


ABSTRACT 

We perform statistical analyses of the infall of dark-matter onto clusters in numeri¬ 
cal simulations within the concordance ACDM model. By studying the infall profile 
around clusters of different mass, we find a linear relation between the mass and 
maximum infall velocities which reach ~ 900 km/s for the most massive groups. The 
maximum infall velocity and the group mass follow a suitable power law fit of the 
form, = (M/mo) 1 - By comparing the measured infall velocity to the linear in¬ 

fall model with an exponential cutoff introduced by Croft et al., we find that the best 
agreement is obtained for a critical overdensity S c = 45. We study the dependence of 
the direction of infall with respect to the cluster centres, and find that in the case of 
massive groups, the maximum alignment occurs at scales r ~ 6/i _1 Mpc. We obtain a 
logarithmic power-law relation between the average infall angle and the group mass. 
We also study the dependence of the results on the local dark-matter density, finding 
a remarkable difference in the dynamical behaviour of low- and high-density particles. 
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1 INTRODUCTION 

The large-scale velocity field provides a significant source of 
information on the distribution of mass fluctuations in the 
Universe. Moreover, the relation between peculiar velocities 
and mass is simple, providing a direct indication of the dis¬ 
tribution of mass independent of assumptions about the bias 
factor of different galaxy types. 

Large scales flows have been addressed both theoret¬ 
ically (Regos & Geller 1989, Bothun et al. 1990, Hiotelis 
2001) and observationally (Dekel 1999, Cecarrelli et al. 2005, 
for a review see Strauss & Willick 1995). On the other hand, 
the dynamical behavior of objects near regions of high den¬ 
sity contrast can be described by the spherical infall model 
(Regos & Geller 1989) with a collapsing streaming motion 
whose amplitude depends only on the distance to the lo¬ 
cal density maximum (Diaferio & Geller 1997). The predic¬ 
tions from the spherical infall model are often used for clus¬ 
ter mass estimates ( e.g. Diaferio 1999). Gunn (1978) and 
Peebles (1976, 1980) derived a linear approximation to the 
infall velocity induced by an isotropic mass concentration 
described by 5(r) = p(r)/pb — 1, where p(r) is the average 
density inside a radius r and pb is the background density. 
The infall velocity Vi„f is then given by 

VlT S = -iffono 6 r$(r), (1) 

where Ho is the Hubble constant and Uq is the cosmolog¬ 


ical density parameter. Yahil (1985) provides a non-linear 
approximation to the exact solution of the infall velocity: 
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however, this equation is not accompanied by a published 
derivation. Regos & Geller (1989), assume that clusters are 
spherically symmetric and obtain a solution for the non¬ 
linear infall velocity of the system as an expansion series, 
which is accurate before orbit crossing takes place and for 
S < 1 — 2. They also quote that this solution is inad¬ 
equate for large density enhancements. Therefore, neither 
their calculations nor Eqs. Qor|2] are expected to hold in 
or near cluster centres, or more generally in virialised re¬ 
gions. In order to overcome this problem, Croft, Dalton & 
Efstathiou (1999) -hereafter, CDE- choose to truncate the 
expression for with an exponential cutoff g-'HD/A 

with S c = 50, which roughly approximates the effect of the 
rapid decrease of infall velocities at small r due to virialised 
motions (see figure 7 in CDE). When adding the exponential 
cutoff, the expression for the linear infall velocity given by 
Eq. Q becomes: 


V/l n f = Hotter S(r)e 5{r)/s *-, 

for the non-linear infall velocity (Eq.|5J this is: 
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CDE claim to obtain a better fit to their simulations using 
this equation, but given the arbitrary nature of this trunca¬ 
tion, they still restrict their analysis to the outer regions of 
haloes (r > 2.5/V 1 Mpc). 

In this paper, we explore the characteristics of stream¬ 
ing motions in the regions surrounding clusters using large, 
high resolution VIRGO N-body simulations, and compare 
these results to predictions from linear theory and the spher¬ 
ical infall model. We will also study the alignments between 
the velocity field and the positions of haloes in the numerical 
simulation, to determine the extent to which a radial approx¬ 
imation to the velocity field around high-density peaks is 
valid. In a related subject, we will provide a further quantifi¬ 
cation of velocity alignments by measuring the fractions of 
particles flowing towards haloes with respect to those stay¬ 
ing at roughly fixed gravitational potentials, and we will 
compare it to the fraction of particles escaping the potential 
wells of haloes. In all these analyses we will also study the 
impact of local density on the infall and alignment of motion 
towards clusters. 

This paper is organised as follows: Section 2 contains 
a brief description of the numerical simulations used in this 
work; Section 3 describes the statistical analyses and the 
resulting dependencies of the infall pattern on cluster mass 
and local density. Finally, Section 4 contains a summary of 
the conclusions drawn from this work. 
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Figure 1. Infall pattern in a slice, 60h~ 1 Mpc a side and I 0h~ 1 
Mpc thick, of the VLS simulation box. The particle density is plot¬ 
ted in a logarithmic gray-scale, smoothed using a top-hat function 
(R = lh~ 1 Mpc). The black solid arrows show the peculiar ve¬ 
locity field. 


2 NUMERICAL SIMULATIONS 

The numerical simulation used in this work was carried out 
by the Virgo Consortium. In particular, we use the pub¬ 
licly available Very Large Simulations, VLS (Yoshida et al. 
2001), characterised by a AC DAI model, 512 3 particles in 
a cosmological box of 480/U 1 Mpc on a side, and a mass 
resolution of 6.86 x 10 w h~ 1 M@. The cosmology is charac¬ 
terised by a matter density parameter S2o = 0.3, cosmolog¬ 
ical constant Ha = 0.7 and expansion rate at the present 
time Ho = 100/t Mpc km/s where h = 0.7. The Cold Dark 
Matter (CDM) initial power spectrum was computed using 
the CMBFAST software (Seljak & Zaldarriaga 1996), and 
is normalised to the abundance of galaxy clusters at z = 0 so 
that ag = 0.9. Dark-matter haloes are identified in the sim¬ 
ulation using a standard friends-of-friends algorithm (Davis 
et al. 1985) using a linking length l, given by l = 0.17n -1 / 3 , 
where n is the mean number density of particles; this link¬ 
ing length corresponds to an overdensity Sp/p ~ 200, which 
characterises the virialized halo and separates it from the 
sorrounding infall material (Cole & Lacey 1995). This large 
simulation contains a large number of halos spanning a wide 
range of masses. The criteria followed in the identification 
of dark-matter haloes leads to ~ 720, 000 identified systems 
with masses ranging from 6.8 x 10 n /i _1 Mg (10 particles) to 
2.35 x 10 15 /i _1 Mq (more than 20, 000 particles). 

In FigureQwe plot a slice of the VLS simulations where 
the distribution of particles is shown in a logarithmic gray¬ 
scale density plot. The dark matter density was smoothed 
using a R = 1 h~ x Mpc top-hat filter and corresponds to a 
slice of 60A -1 Mpc a side and 10fi _1 Mpc deep. The black 
arrows indicate the projected velocity field of 20% of the 


particles in the slice. It can be clearly appreciated the sys¬ 
tematic infall flow onto dense regions, as well as outflow 
patterns from the lower density regions. 

We calculate the halo peculiar velocity by averaging the 
individual group member velocities. In Figure |2] we show 
the normalised distribution function of the peculiar veloc¬ 
ity of halos with masses between 10 11 h _1 Mq — 10 12 fi _1 Mg 
(dotted line), 10 12 /i _1 Mg — 10 13 /i _1 Mq (short dashed line), 
10 13 /i _1 Mq — 10 14 /i _1 Mq (long dashed line) and more 
than 10 14 /i _1 Mq (solid line); and the distribution func¬ 
tion of peculiar velocities of the dark matter particles (thick 
short-dashed line). As can be seen, even the most massive 
groups have similar mean-velocities than dark-matter parti¬ 
cles (Einasto et al. 2005, Shaw et al. 2005). However, only 
dark-matter particles show a high velocity tail (V > 1500 
km/s), whereas groups rarely reach peculiar velocity values 
above V = 1000 km/s (Shaw et al. 2005). 


3 VELOCITY FIELD ANALYSIS 

In this section we study the statistical properties of the ve¬ 
locity field in the neighborhoods of dark-matter haloes, up 
to scales of 40/V 1 Mpc. We begin by analysing the mean 
infall velocity of particles onto groups, and the distribution 
of the peculiar velocity orientation of the particles in the 
outskirts of the groups (Subsection 13.111 . In Subsection IQ 
we study the dependence of the results on the local dark- 
matter density. We perform our analysis for four diferent 
group mass samples (see Table 0. 

All our analises are performed using individual parti¬ 
cles in the simulation as opposed to the recent use of sub¬ 
structure identified in simulations in several works (see for 
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Figure 2. Normalised distribution function, /„, of the velocity 
module V = \J v%. + for particles (thick short dashed line) 

and groups with masses in the ranges 10 n /! _1 Mg — 1 0 1 2 h ~ 1 A/q 
(dotted line), 10 12 /i — 1 Mq — 10 13 /i -1 Mq (short dashed line), 
10 13 /i _1 Mq — 10 14 /i _1 Mg (long dashed) and more than 10 14 
h~ 4 Mq (solid line). 


Figure 3. Average infall velocity onto groups in samples SI, S2, 
S3 and S4 (dot dashed, short dashed, long dashed and solid lines, 
respectively) as function of scale. Dotted lines represent the re¬ 
sults from the linear fit with an exponential cutoff (Eq. [o] with 
8 C = 50. 


Table 1. Description of our samples of dark-matter haloes cor¬ 
responding to different ranges of mass and the number of groups 
identified. 


Sample 

Mass interval [h 1 Mq] 

Number of groups 

SI 

3.4 x 10 12 - 6.8 x 10 12 

73509 

S2 

6.8 x 10 12 - 6.8 x 10 13 

72465 

S3 

6.8 x 10 13 - 6.8 x 10 14 

6097 

S4 

> 6.8 x 10 14 

92 


instance Benson, 2005, and Wang et al., 2005). We consider 
that our studies should not change importantly if we were 
to repeat the following analyses using substructure since we 
are particularly interested in regions outside the halo cen¬ 
tres. At such distances, we expect the average velocities of 
individual particles and substructre to be similar. This is 
opposed to the clear differences inside dark-matter haloes, 
where dynamical friction would slow down the center of mass 
of a sub-halo. 


3.1 The average velocity infall 

We start our analysis of the velocity field by computing the 
group-averaged infall velocity as follows: we select groups in 
a specific mass range, and calculate the mean particle infall 
velocity as a function of distance to the group centres. In 
Figure |3] we plot the mean infall velocity profile for the four 
subsamples SI, S2, S3 and S4 (dotted, short-dashed, long- 
dashed and solid line respectively) described in Tabled As 
it can be seen, the maximum infall velocity, V£tf X , increases 


with the group mass, with values V^f X ~ 250 km/s for 
groups in sample SI, and V^f x = 950 km/s for those in S4, 
the most massive groups in our analysis. The rapid decrease 
of infall velocities at small scales indicates the boundary of 
the virialised region. It can also be seen that the scale of the 
maximum infall velocity correlates well with group mass, 
a somewhat expected effect given the correlation between 
mass and virial radius. We notice a more rapid decrease of 
the infall pattern for the most massive groups compared to 
the gentle decline shown by low mass systems. We compare 
the actual infall velocity measured in the simulations with 
the prediction from the CDE model. In order to do this, we 
calculate the averaged density enhancement inside a sphere 
of radius r centered in each group and derive the linear in¬ 
fall velocity using an exponential cutoff given by Eq. [3] for 
the same mass intervals described in Table Q The results 
are plotted in Figure [2] where it can be appreciated that 
although this model provides a reasonable agreement to the 
infall around high mass haloes, there are still large discrep¬ 
ancies at small and large scales, specially for the lower mass 
ranges explored here. 

The results shown in Figure ^indicate a tight relation 
between V^f x and group mass in both, the patterns mea¬ 
sured in the simulations, and the results from Eq.£3 In order 
to address the infall correlation with group mass we mea¬ 
sure the mean infall profile for different mass ranges, and 
compute V™f x ■ In Figure El we show the maximum mean 
infall velocity as a function of group mass in the simulation 
(solid circles), and from Eq.|3](open triangles). 

A suitable power-law fit to the simulation results is: 


V -max 
inf 
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Figure 4. Maximum infall velocity as a function of group mass 
measured in the simulation (black dots) and its corresponding 
power-law fit (solid line). The open triangles show the same re¬ 
lation as obtained from the maxima in Eq. |3] We also show for 
comparison the results from Eq. [7| assuming different values of 
S c = 40, 45, 50, 55, 60 (dotted lines). The best-fit to the simulation 
results is given by the dashed line, which corresponds to <5 C = 45. 


where the parameter values are mo = 1.29 x 10 5 , and 7 = 
0.29. Using the linear fit with the exponential cutoff given 
by Eq. m it is possible to infer a relation between the scale 
of the maximum infall velocity and the halo mass: 


max 
T inf 



1/3 

M 1/3 . 


( 6 ) 


where the factor A is defined as function of 5 C as A = 
V(2 + 3/<5 C ) 2 + 12/8c - 2 - 3/6 C , with A > 0 for any given 
positive value of <5 C . Replacing r by r" T l , a j x in Eq. [3] we find 
the expression for V/£f X : 


Vinf (r ri 


(9tt) 


- 1/3 




2 -A 

~AP3~ 


e^M 1/3 . (7) 


By setting Ho = 70, fi m = 0.3 and <5 C = 50, we obtain 
Vi n f(rm ax ) ~ M 1 ^ 3 / 2.25 x 10 4 , which is very similar to our 
power-law fit in Eq.|5] We now calculate V/£f X (M) from Eq. 
Qusing different values of 5 C = 40, 45, 50, 55 and 60, and plot 
the results in Figure2](dotted lines). As it can be seen in this 
figure, the best fit to the simulations is achieved when using 
5 C ~ 45 (dashed line). The power-law behaviour of 
brakes down at low mass halos, when M < 10 13 ft _ 1 Mg. At 
this point, the maximum infall velocity does not change con¬ 
siderably with halo mass, as can be seen in the infall pat¬ 
terns for samples SI and S2 in figure 03 The reason behind 
this may be related to the flattening of the bias between 
dark-matter haloes and mass b(M), which occurs at similar 
mass scales; namely, the clustering around groups of masses 
M < 10 13 /i _1 Mq can be considered to be almost indepen¬ 
dent of halo mass. 

In order to characterise the alignments of the infall of 


Figure 5. Average cosine of the angle between the cluster-centric 
position and peculiar velocity of dark-matter particles around 
haloes in the simulation. Groups are divided into samples span¬ 
ning different mass ranges indicated in Table [T] (samples SI, S2, 
S3 and S4, in dotted, short dashed, long dashed and solid lines, 
respectively). The inset shows the relation between {eos{6)) max 
and group mass. The dotted line corresponds the best log-linear 
fit to this relation, cos{9) rna x = 'ylog(MfCo), where 7 = 0.17 and 
Co = 5.18 x lO 1 !). 


dark-matter in the outskirts of haloes, we centre our origin 
on each halo, and measure the angle 9 between the particle 
position vector and its peculiar velocity. In order to detect 
variations with the halo mass, we analyse subsamples SI, S2, 
S3 and S4 described in Table Q Figure 03 shows the results 
from averaging ( cos(9 )) over each sample of haloes (dotted, 
short dashed, long dashed and solid lines for samples SI, S2, 
S3 and S4, respectively). 

Dark-matter particles infalling towards S4 groups show 
the maximum alignment, cos{9) rnax ~ —0.77 at r ~ 5.5/i _1 
Mpc. In contrast, the peculiar velocity of particles in the 
outskirts of groups in sample SI show the lowest alignment, 
with cos{6) max ~ —0.3 at r ~ l/i -1 Mpc. As can be seen, 
we obtain a maximum ‘coherence’ for the alignment in a rel¬ 
atively small range of scales: 1 ft -1 Mpc < r < 6 / 1 -1 Mpc. 
The velocity field of the particles around the most massive 
groups (sample S4) show a significant alignment at large 
scales (r > 30h -1 Mpc). However, for less massive systems, 
( cos{9 )) is significantly reduced. It can also be appreciated 
that the maximum alignment between the cluster-centric po¬ 
sition and the peculiar velocity of the particles is a function 
of group mass. We show this relation in the inset of Figure 
13 for different masses (same samples as those in Figure 0. 
An empirical logarithmic power-law fit, 

cos(9)max(M) =-ylog{M/C 0 ), (8) 

with parameters 7 = 0.17 and Co = 5.18 x 10 10 , gives a very 
good agreement to the simulation results, as can be seen in 
dotted line in the inset of Figure 0 

We perform a further study of the infall-angle in the 
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outskirts of dark-matter haloes by comparing the number 
of particles in each shell of radius r with peculiar veloci¬ 
ties pointing outwards from the group center, to the num¬ 
ber pointing towards the group center. We measure Ni(r), 
N m {r) and N 0 (r), which indicate the number of particles 
characterised by 1/3 < cos(9) (infall), —1/3 < cos(9) < 1/3 
(intermediate) and cos(9) < — 1 /3 (outflow), respectively, at 
a distance r from the group centre. We average across our 
sample, and calculate the ratios Di(r ) =< ( Ni(r)/N m (r )) > 
and D 0 (r) =< ( N 0 (r)/Nm[r )). Figure[S]show the results for 
Di(r) (thick lines) and D 0 (r ) (thin lines) for the samples SI, 
S2, S3 and S4 (dotted, short-dashed, long-dashed and solid 
lines, respectively). As can be seen, samples SI and S2 show 
D 0 > 1 at r > 30/i _1 Mpc, indicating that the large density 
inhomogeneities marked by the halo centres are still corre¬ 
lated with the distribution of mass at such large distances. 

The fact that the maximum of Di and minimum of D 0 
coincide with one another, as well as with the maximum 
alignment as measured by cos(9), is to be expected since the 
maximum alignment occurs when the fraction of infalling 
particles reaches its maximum (compare figure 0 and |SJ|. 
We also notice that these maxima occur at larger distances 
for higher mass haloes, consistent with the maximum in¬ 
fall velocities occurring further away from the most massive 
haloes. We notice that the maximum infall velocity occurs 
~ 30% closer to the group than the maximum velocity align¬ 
ment. This suggests that alignments are a better choice than 
infall velocities for the determination of the distance where 
the halo virialisation breaks the trend of increasing infall 
and velocity alignments toward the halo centres. In general, 
our results are consistent with high mass haloes being more 
relaxed systems with isotropic velocity distributions. 

3.2 The infall dependence on local density. 

In this subsection, we analyse the dependence of the infall 
pattern, the infall-angle and the parameters Di and D 0 on 
the local density of the dark-matter particle being consid¬ 
ered in the analysis. We define a local number density, p n , 
as the number of particles in cubic cells of side l c = 3/i -1 
Mpc (note this is not a proper density definition, just a 
simple parameter to characterise the environment of a dark- 
matter particle), and a mean numerical density, p„, as the 
mean number of particles expected in cells of volume V c = lc 
(in this particular numerical simulation, and for the chosen 
value of lc, the mean numerical density is p n = 33). All par¬ 
ticles within a cubical cell are assigned the same density and 
we divide our analysis according to the following ranges of 
local density: high, intermediate and low densities defined as 
Pn/pn > 2, 1/2 < p n /p n < 2 and pn/pn < 1/2, respectively. 

In order to test the possibility of a further dependence of 
the infall pattern on mass and local density, we analysed the 
fraction of particles laying in diferent density ranges as func¬ 
tion of scale. The fraction of particles in high, intermediate 
and low density regions are calculated using fhi = Nhi/Ntot, 
fav = Nav/Ntot and fi a = Ni 0 /N t0 t, where N h i, N av and N io 
are the number of particles in high, intermediate and low 
density regions respectively, and N to t = Nhi + N av + N 0 is 
the total number of particles within each distance interval. 
In figure Q we plot with thick lines the fraction of particles 
fhi (solid lines), f a v (dashed lines) and fi a (dotted lines) for 
the S4 M > 6.8 x 10 14 sample, and with thin lines the corre- 



r [Mpc h _l ] 


Figure 6. The average ratio between the number of particles with 
1/3 < cos(9) < 1 ( N out fi ow ) and -1/3 < cos(9) < 1/3 (Nm), 
referred to as Di (black); and the ratio between the number of 
particles with —1 < cos(9) < —1/3 (.%:„/«//) and N m , refereed 
to as D a (red). These ratios are plotted for the different mass 
ranges used in previous plots: 3.45 X 10 11 — 6.8 X 10 12 /i _1 Mq 
(dotted line), 6.8 X 10 12 — 6.8 X 10 13 /i _1 Mq (short dashed line), 
6.8 X 10 13 — 6.8 X 10 14 /i -1 Mq (long dashed line) and M > 6.8 X 
10 14 /i 1 AIq (solid line). 

sponding fraction of particles for sample SI (M < 6 . 8 x 10 12 ). 
In this figure, we can appreciate that the infall patterns cor¬ 
responding to fhi for the two mass range sample differ in 
less than 20 % in scales r < 20 Mpc/h while the patterns of 
fav and fi 0 are similar within 10 % in the same scales. 

In Figure [H] we plot the average infall velocity of parti¬ 
cles with p/p > 2 (solid lines), p/p ~ 1 (dashed lines) and 
p/p < 1/2 (dotted lines) as a function of the scale, for the 
four group samples SI, S2, S3 and S4 (panels A, B, C and 
D respectively). It can be noticed that more massive haloes 
tend to show less differences in the infall patterns of high 
and low density particles (almost no difference in panel D), 
whereas low mass groups (panel A) show a difference be¬ 
tween maximum infall of high and low density particles of 
~ 200 km/s (more than 60% of the maximum infall for the 
high density particles). We provide fits to the vs. mass 

relation for the different ranges of local density in table [ 5 ] 
where it can be seen that both the slope, 7 , and mo increase 
for decreasing local density values. Note that the global fit 
is similar to that of p n > 2pn, which indicates that the sig¬ 
nal from high density regions dominates the overall infall 
velocity. 

In order to get a deeper insight on the flow patterns 
in the outskirts of groups, we also study the dependence of 
the infall-angle on density; we compute the average ( cos(9 )) 
for particles characterised by a local density p n . Figure [5] 
shows these results for the mass samples defined in Ta¬ 
ble Q (panels A, B, C and D for samples SI, S2, S3 and 
S4 respectively) for particles in high density environments, 
p » p (solid lines); intermediate densities, p ~ p (dashed 
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Figure 7. Fraction of particles in high, intermediate and low den¬ 
sity regions , _/),, (solid lines), f av (dashed lines) and fi 0 (dotted 
lines), respectively as a function of scale. The thick lines corre¬ 
spond to the most massive groups (sample S4) and the thin lines 
to the least massive ones (sample SI). 


Figure 8. Mean infall velocity as a function of scale. Halo sam¬ 
ples correspond to SI, S2, S3 and S4, and are defined in Table 
0 (panels A, B, C and D). The different lines show the average 
infall patterns of particles in regions with density p » p (solid 
line), p ~ p (dashed line) and p << p (dotted line). 


Table 2. Values obtained for the parameters 7 and 7 ?i 0 in the fit Table 3. Values obtained for the parameters 7 and Co in the fit of 

of vs. mass, for different values of local density. cos(6) m ax vs. mass, considering different ranges of local density. 


Local density parameter ( p n ) 

7 

mo 

All pn values 

0.295 

1.29 x 10 s 

pn 2 p-n 

0.281 

4.22 x 10 4 

pn/'Z < rho n < 2 Pn 

0.404 

4.85 x 10 7 

pn < pn/2 

0.442 

1.36 x 10 s 


line) and low densities, p « p (dotted line). Each group 
sample shows diferent features in their alignment profiles 
including different maximum alignments; we notice that for 
higher mass groups the maximum alignment is more signifi¬ 
cant, as can be seen in Figure^] For the least massive groups 
(sample SI in panel A) we obtain a maximum alignment of 
( cos(6 )) ~ —0.5 at r ~ 1 Mpc for low density parti¬ 

cles, and a lower maximum alignment for particles in high 
density regions, ( cos(6)) ~ —0.25 at r ~ 4 h _1 Mpc. How¬ 
ever, at larger scales, r > 6 h~ x Mpc, high-density particles 
are better aligned than low-density particles. Samples S2 
and S3 (intermediate masses) show a similar behaviour: the 
maximum alignment of low-density particles around groups 
in sample S2 reaches ( cos{6 )) ~ —0.6, whereas the high- 
density particles reach ( cos(9 )) ~ —0.4, both at r > 8he 1 
Mpc. As can be seen, the peculiar velocities of higher-density 
particles are more aligned with the direction to the centre of 
the halo. In the case of group sample S3 the particles in low 
density regions have a maximum alignment ( cos(9 )) ~ —0.8 
whereas high-density particles show a maximum alignment 
cos{6 ) ~ —0.6. Again, at scales r > 20hA 1 Mpc, higher den¬ 
sity particles are more aligned. In group sample S4 it can 


Local density (pn) 

7 

Co 

All p n values 

0.176 

5.18 x 10 10 

Pn > 2 pn 

0.178 

6.16 x 10 10 

Pn/ 2 < rho n < 2 p n 

0.159 

2.66 x 10 9 

pn < Pn / 2 

0.157 

1.05 x 10 9 


be observed that particles in low and high density environ¬ 
ments show the same maximum alignment cos(d) ~ —0.8 
at r ~ 6h _1 Mpc, while particles with intermediate densi¬ 
ties are more strongly aligned, with cos(6) ~ —0.9. For this 
particular group sample, we do not find a significant dif¬ 
ference in the alignment of high, medium, and low density 
particles. Table [3] shows the parameter values for the fits to 
cos(9)max vs. mass for the different ranges of local densities 
and halo mass considered here. As can be seen, both the 
slope 7 and Co are larger for higher densities. This indicates 
that particles in denser regions are moving in a direction 
better aligned to the center of the groups. As in the case 
of the maximum infall velocity, it can also be seen that the 
parameter fits resulting from all the particles surrounding 
groups is similar to the fit from the particles in the densest 
regions, indicating that the latter dominate the alignment 
signal. 

We now search for possible dependencies of the frac¬ 
tion of infalling and outflowing particles, D 0 and Di defined 
above, on the local density of dark-matter particles. Figure 
HOl presents these results for our four group samples (panels 
A, B, C and D for samples SI, S2, S3 and S4, respectively) 
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Figure 9. Dependence of { cos(6 )) on local density for our four 
mass groups samples. Panel A: ( cos(8 )) for groups in mass sample 
SI, computed for particles in regions with p >> p (solid line), 
p ~ p (dashed line) and p « p (dotted line). Panel B: same as 
A, for groups in sample S2. Panel C: r groups in sample S3. Panel 
D: groups in sample S4. 


for particles in low- (dotted line), intermediate- (dashed line) 
and high-density (solid line) regions. As can be seen, the ra¬ 
tio Di is larger for low-density particles at small separations, 
whereas at larger scales, the trend is similar for particles in 
both, high and low density regions. Inspection of panel D, 
indicates that an important difference of three orders of mag¬ 
nitude can be achieved between the number of infalling to 
outflowing particles for high density regions. Note that this 
maximum occurs at a scale ~ 50% larger than that in low 
density regions. In high density regions, however, this dif¬ 
ference is not so important. It should also be noticed that, 
although not detected in Figure 6 for the overall ratios of 
infalling and outflowing particles, particles in the lower den¬ 
sity regions show a minimum in Di at slightly (marginally) 
larger distances from the halo centres than the correspond¬ 
ing maximum in D 0 - 


4 CONCLUSIONS 

We now summarise the main conclusions we arrived at dur¬ 
ing the series of analyses performed on the VLS simulations. 
We studied the peculiar velocity field around dark-matter 
haloes of varying masses, considering infall velocities, the 
infall angle with respect to the halo centres, and the frac¬ 
tions of infalling and outflowing particles. We also studied 
variations in these quantities as the local dark-matter par¬ 
ticle density varies. 

We now list our main results and conclusions: 

• The maximum infall velocity around the most mas¬ 
sive groups identified in the numerical simulation (~ 
10 15 /i _1 Mq), is Vinf ~ 900 km/s. The results for the lowest 



Figure 10. Dependence of N ou tflow/Nm (black) and 
Ninfall/Nm (gray) with the local density of dark-matter 
particles, for four subsamples of halos (panels A, B, C and D for 
samples SI, S2, S3 and S4, respectively). Results for Dj and D 0 
for particles characterised by p/p > 2 are shown in solid lines, 
whereas dashed lines shows the results for p ~ p, and dotted 
lines the results for p/p < 1/2. 


mass group sample (~ 10 12 h^ 1 Mq) indicate a maximum in¬ 
fall velocity of ~ 220 km/s, decreasing very slowly at larger 
distances from the halo centres. 

• The maximum infall velocity occurs in a narrow range 
of scales, 2/U 1 Mpc < r m ax < 6h _1 Mpc even when con¬ 
sidering a broad range of group masses (10 12 /i _1 Mq to 
10 15 /i -1 M q ). 

• The exponential cutoff linear model analysed by Croft 
et ah, works better for the most massive groups at small 
scales. It does not reproduce the infall pattern for groups 
samples with M < 10 14 Ii~ 1 Mq, neither the large scale infall 
signature onto massive systems in any of our subsamples SI 
to S4. 

• We find a power law relation between the maximum 
infall velocity and mass. The Croft et al. model with an 
exponential overdensity cutoff provides a good fit for <5 C ~ 
45. 

• Large mass groups show a stronger alignment between 
the dark-matter particle peculiar velocity and its cluster¬ 
centric position than do low mass groups. Regardless of the 
halo mass, this maximum alignment occurs at intermediate 
scales, 2h~ 1 Mpc < r m ax < 8/i _1 Mpc. 

• There are remarkable differences between the dynam¬ 
ical behaviour of low- and high-density particles. In gen¬ 
eral, higher density particles exhibit higher infall velocities 
and lower alignment than particles in lower density environ¬ 
ments. The difference in infall is less evident when consider¬ 
ing high mass haloes, and the alignment becomes similar for 
high and low density particles for larger distances to halo 
centres. However, at large enough distances, high-density 
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particles become even more aligned than lower density par¬ 
ticles, specially for low mass haloes. 

We have also considered the possibility that substruc¬ 
ture could affect our analysis. According to the results of 
Aubert, Pichon & Colombi (2004) the infall pattern of the 
largest halos is affected by substructure only at small scales 
(r < 2 Mpc). Since the maxima occurs well beyond this scale, 
we are confident that our results are not likely to change un¬ 
der the presence of substructure. 

In spite of the fact that mass flows toward clusters in 
an anisotropic fashion preferentially through the filamentary 
structure, the results provided give useful mean values of ve¬ 
locity flows which can be used to asses the systematic infall 
onto growing structures. The particular analysis performed 
provide direct measures of the maximum infall pattern and 
halo mass dependence and can be used to estimate gas dy¬ 
namical processes such as ram pressure stripping that may 
initiate galaxy transformation in the outskirts of clusters. 


ACKNOWLEDGMENTS 

MP acknowledges receipt of a CONICET fellowship. NDP 
was supported by a Proyecto FONDECYT Postdoctoral 
3040031. This work was supported in part by the ”ESO 
Center for Astrophysics” at Catolica and by the European 
Commission’s ALFA-II programme through its funding of 
the Latin-american European Network for Astrophysics and 
Cosmology, LENAC. The simulations in this paper were 
carried out by the Virgo Supercomputing Consortium us¬ 
ing computers based at the Computing Centre of the Max- 
Planck Society in Garching and at the Edinburgh Paral¬ 
lel Computing Centre. The data are publicly available at 
http://www.mpa-garching.mpg.de/NumCos This research 
has made use of NASA’s Astrophysics Data System. 

REFERENCES 

Aubert D., Pichon C. & Colombi S., 2004, MNRAS, 352, 344. 
Bothun G.D., West M.J., Mould J.R. & Schommer R.A., 1990, 
ApJ, 353, 344. 

Ceccarelli M.L., Valotto C., Padilla N., Lambas D.G., Giovanclli 
R. & Haynes M., 2005, ApJ, 622, 853. 

Cole S. & Lacey C., 1996, MNRAS, 281, 716. 

Croft R.A.C., Dalton G.B. & Efstathiou G., 1999, MNRAS, 305, 
547. 

Davis M., Efstathiou G., Frenk C.S. & White S.D.M., 1985, ApJ, 
292, 371. 

Dekel A., 1999, Formation of Structure in the Universe, ed. Dekel 
A. and Ostriker J.P. (Cambridge University Press), p. 250. 
Diaferio A., 1999, MNRAS, 309, 610. 

Diaferio A. & Geller M.J., 1997, ApJ, 481, 633. 

Einasto M., Suhhonenko I., Heinamaki P., Einasto J., & Saar E., 
2005, A&A, 436, 17. 

Gunn J.E., 1978, Observational Cosmology, 8th Annual Saas-Fee 
Course, edited by A. Maeder, L. Martinet and G. Tammann 
(Geneva Observatory, Sauverny). 

Hiotelis N., 2001, A&A, 375, 338. 

Peebles P.J.E., 1976, ApJ, 205, 318. 

Peebles P.J.E., 1980, Large Scale Structure of the Universe 
(Princeton University, Princeton). 

Regos E. & Geller M.J., 1989, ApJ, 98, 755. 

Seljak U. & Zaldarriaga M., 1996, ApJ, 469, 437. 


Shaw L., Weller J., Ostriker J.P. & Bode P., astro-ph/0509856 
Strauss M.A. & Willick J.A., 1995, Phys. Rep., 261, 271. 

Yahil A.. 1985, The Virgo Cluster, ed. Richter O. and Binggeli B. 
(ESO, Garching), p. 359. 

Yoshida N., Sheth R., & Diaferio A., 2001, MNRAS, 328, 669. 


